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ABSTRACT 

Adopting  appropriate  postures  during  manual  material¬ 
handling  tasks  is  the  key  to  reducing  human  joint  injuries. 
Although  much  experimentation  has  been  conducted  in  an  effort 
to  model  lifting,  such  an  approach  is  not  general  enough  to 
consider  all  potential  scenarios  in  material  handling.  Thus,  in 
this  paper  an  optimization-based  motion  prediction  method  is 
used  to  simulate  realistic  lifting  postures  and  predict  joint 
torques  to  evaluate  the  risk  level  of  injury. 

A  kinematically  realistic  digital  human  model  has  been 
developed  such  that  the  complicated  musculoskeletal  human 
structure  is  modeled  as  a  combination  of  serial  chains  using  the 
generalized  coordinates.  Lagrange’s  equations  of  motion  and 
metabolic  energy  rate  are  derived  for  the  digital  human. 

The  proposed  method  has  been  implemented  to  predict  and 
evaluate  the  lifting  postures  based  on  the  metabolic  rate  and 
joint  torques.  Our  results  show  that  different  amount  of  external 
loads  and  tasks  lead  to  different  human  postures  and  joint 
torque  distribution,  thus  different  risk  level  of  injury. 

Keywords:  lifting,  posture,  motion  prediction,  joint  torque, 
Lagrangian,  optimization 

INTRODUCTION 

Digital  human  models  allow  extensive  analysis  and 
simulation  to  be  carried  out  in  a  virtual  world  rather  than  with 
an  expensive  mock-up.  The  users  are  able  to  get  feedback  from 
the  simulation  of  certain  tasks  and  performances  of  a  human. 
For  material  handling,  the  simulation  feedback  is  an  important 
aspect  of  the  task  design  procedure.  The  goal  is  to  create 
methods  and  tools  that  allow  for  human  analysis,  design,  and 
evaluation,  without  actually  having  a  physical  mock-up  of  the 
environment. 


This  paper  addresses  one  such  tool  that  will  allow  a  user  to 
predict  human  motion  and  performance  level  for  a  given  task. 
To  that  end,  we  propose  a  general  approach  to  mathematics- 
based  dynamic  motion  prediction  and  analysis.  We  assume  that 
a  human  naturally  moves  in  such  a  way  as  to  minimize  certain 
cost  functions.  In  the  case  of  human  behavior,  we  call  the  cost 
functions  the  human  performance  measures.  Therefore,  we 
propose  the  implementation  of  an  optimization-based  approach 
whereby  human  performance  measures  are  developed  and 
utilized  to  predict  the  realistic  motion. 

During  manual  lifting  tasks,  the  lumbar  spinal  column 
(low-back)  is  the  region  where  most  back  injuries  occur  as  it  is 
subjected  to  much  greater  loads  than  any  other  area  of  the  spine 
(Gallagher  et  al  2002).  The  reason  for  the  high  forces 
produced  on  the  lower  back  during  lifting  activities  is  the 
vigorous  contraction  of  back  muscles  to  counteract  the  weight 
being  lifted.  Gallagher  et  al  have  experimentally  estimated  the 
low-back  muscle  forces  for  four  different  cable-lifting  postures 
(kneeling  on  one  knee,  kneeling  on  two  knees,  stooping,  and 
standing)  from  the  EMG  data.  The  lumbar  moments  for  each 
case  were  calculated  from  the  muscle  forces  and  muscle 
moment-arms.  These  were  also  compared  with  the  cases  for 
increased  cable  load.  Chaffin  (1997)  has  developed  two-  and 
three-dimensional  computerized  human  models  for  simulation 
of  static  strength  during  work.  This  method  can  also  predict 
some  critical  values  of  low-back,  such  as  back  compression 
loads  and  disc  shear  loads. 

Considerable  research  has  been  done  in  postural  analysis  in 
material  handling  (Bobick  1987;  Nussbaum  and  Torres  2001; 
Gallagher  1987;  Gallagher  et  al  2002;  Chung  and  Kee  2000). 
Most  of  them  used  experiment-based  methods. 
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On  the  other  hand,  there  are  several  studies  on  human 
posture  and  motion  simulation  based  on  optimization.  One  of 
the  most  widely  used  cost  functions  in  human  simulation  studies 
is  the  energy  expenditure.  Hase  and  Yamazaki  (1997)  have 
simulated  walking  and  rowing  motions  using  rigorous  muscle 
energy  expenditure  formulas  and  human  model.  They  have 
used  the  Newton-Euler  method  to  calculate  the  joint  torques. 
However,  calculating  the  muscle  energy  expenditure  directly 
from  the  element  muscle  fiber  model  requires  tremendous 
computational  efforts.  Umberger  et  al  (2003)  formulated 
rather  simplified  yet  exact  muscle  energy  expenditure,  based  on 
Hill-type  muscle  model  and  coefficients.  Some  simple 
examples,  such  as  isolated  muscle  actions,  single  joint  motion, 
and  locomotion,  were  simulated  and  the  computed  energy 
values  were  compared  with  experimental  results.  Another 
approach  for  obtaining  energy  cost  has  been  used  by  Alexander 
(1997).  Alexander’s  energy  formula  is  very  simple  and 
experiment-based.  That  energy  cost  was  minimized  to  solve  for 
realistic  human  arm  trajectories. 

The  energy  expenditure  can  also  be  related  to  fatigue 
prediction.  Fatigue  is  a  complex  phenomenon,  and  it  may  be 
determined  by  several  factors,  including  psychological  health  of 
the  subjects.  Although  much  is  known  about  fatigue,  its  exact 
cause  remains  to  be  determined.  At  the  biomechanical  level, 
many  literatures  explain  fatigue  as  loss  of  energy  (Nicolson 
2003).  In  the  human  gait  motion,  Anderson  and  Pandy  (2001) 
suggested  that  minimizing  muscle  fatigue  at  each  instant  is 
roughly  the  same  as  minimizing  metabolic  energy  expended  per 
unit  distance  traveled  over  the  duration  of  the  gait  cycle.  In 
fact,  it  is  well  known  that  energy  expenditure  and  muscle 
fatigue  have  positive  correlation  (Sahlin  et  al  1998,  Khang  and 
Zajac  1989).  Therefore,  minimum  metabolic  energy 
expenditure  indicates  less  muscle  fatigue. 

In  this  paper,  we  first  present  a  modeling  method  for 
representing  human  upper  body  that  can  also  be  used  for 
modeling  any  kinematic  chain.  Muscle  models  are  also 
considered  for  dynamics  formulation.  We  then  derive  general 
equations  of  motion  and  human  energy  expenditure  formula  in 
joint  space.  Using  these  equations,  an  optimization  formulation 
is  developed  that  allows  for  a  computer  program  to  simulate  the 
natural  human  motions,  joint  torques  and  energy  level.  Then, 
some  examples  of  motion/posture  prediction  for  manual  lifting 
tasks  followed  by  joint  torque  evaluations  are  illustrated.  By 
predicting  the  joint  torques  at  the  back  and  shoulder,  the  risk 
level  of  injury  can  be  evaluated.  Conclusions  and  future  work 
are  addressed. 

KINEMATIC  HUMAN  MODELING  AND  GENERALIZED 
COORDINATES 

To  describe  the  human  postures  and  motions,  we  use  a 
kinematic  pair  as  used  in  the  study  of  robotics.  A  digital  human 
model  called  Santos™,  developed  by  The  University  of  Iowa, 
has  more  than  100  degree-of-freedom  (DOF),  where  all  the 
joints  are  modeled  as  re  volute  joints.  The  shoulder  joint,  for 
example,  has  three  degrees  of  freedom  and  is  usually  modeled 


as  a  spherical  joint  in  the  literatures.  In  our  model,  we  combine 
three  re  volute  joints  at  one  location  to  represent  the  three 
degrees-of-freedom  spherical  shoulder  joint.  For  upper  body 
motion  prediction,  we  use  a  21  DOF  Santos™  model  from  the 
waist  to  the  right  hand  as  shown  in  Figure  1.  Here,  joint  1 
corresponds  to  lateral  bending  of  low-back,  joint  2  corresponds 
to  flexion/extension  of  low-back,  and  so  on. 


The  position  vector  of  a  point  of  interest  on  the  end- 
effector  of  a  human  articulated  model  can  be  written  in  terms  of 
joint  coordinates  as 

x  =  x(q)  (1) 

where  q  =  [qx  ...  eR"  is  the  vector  of  ^-generalized 

joint  coordinates  defining  the  motion  of  a  link  with  respect  to 

another.  The  global  position  vector  x(q)  can  be  obtained  from 

the  multiplication  of  the  4x4  homogeneous  transformation 

matrices  7_1T7 ,  defined  by  the  Denavit-Hartenberg  (D-H) 

representation  method  (Denavit  and  Hartenberg,  1955). 

cosg,  -cos  a,  sing,  sin  or,  sing,  ajcosqi 

,  sing,  cos  a,  cos  g.  -sin  <x  cos  g,  a.  sing,. 

^.(g,)  =  Hl  1  Hl  1  1  Hl  (2) 

0  sin  al  cos  at  dl 

0  0  0  1 


where  g.  is  the  joint  angle  from  X7l  axis  to  the  xz  axis  for  a 
re  volute  joint,  dt  is  the  shortest  distance  between  xj_l  and  x7 
axes,  af  is  the  offset  distance  between  z7.  and  zi_l  axes,  and 
a 7.  is  the  offset  angle  from  Z7  l  and  Z7  axes.  Then  the  4x4 
transformation  matrix  °T7  used  to  represent /th  joint  coordinate 
system  with  respect  to  the  global  base  coordinate  system  (0th)  is 
0 1) (q)  =°  T,  (q, )‘ T2 (q2 T;. (q, )  (3) 

We  use  the  augmented  4x1  vector  °r.  and  7r.  to  express  the 

Cartesian  coordinate  of  the  point  fixed  in  the  /th  local  frame  in 
terms  of  global  and  local  coordinate  system,  respectively: 


'x(q)' 

r  - 

~zx~ 

1 

?  */■ 

_  1 

where  l\  is  the  fixed  point  at  the  link  i  and  expressed  with 
respect  to  the  zth  coordinate  system.  Using  these  relationships, 
°r.  can  be  written  as 

°i*  =°Tl(ql,...,qiyri  (5) 
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MUSCLE  MODEL  AND  JOINT  STIFFNESS 

A  muscle  is  perceived  as  an  actuator  input  into  the 
musculoskeletal  system  where  each  joint  is  powered  by  one  or 
more  such  actuators.  In  other  words,  muscle-generated  forces 
empower  and  control  all  human  movements  and  postures.  The 
muscle-tendon  systems  have  certain  elastic  properties  that 
generate  restoring  forces  at  the  joint.  These  restoring  forces 
need  to  be  included  in  the  force  equilibrium  of  muscle-link 
system. 

Several  muscle  models  were  proposed  in  the  literature, 
perhaps  most  notably  the  Hill’s  model  (Hill,  1938)  and  the 
Zajac  Muscle  model  (Zajac,  1989),  which  is  a  modified  Hill’s 
model.  Most  of  the  proposed  muscle  models  have  two  main 
components:  contractile  components  and  series/parallel  elastic 
components  (Figure  2). 


series  elastic  contractile 
component  component 


F 


parallel  elastic  component 


F 


Figure  2.  A  typical  Muscle  Model 


The  muscle  contractile  elements  generate  tension  force  by 
contracting,  and  act  as  actuators.  The  elastic  elements  of 
muscles  generate  restoring  forces  in  the  corresponding 
independent  degree-of-freedom  for  each  joint  motion  (where 
each  joint  may  have  many  muscle  degrees  of  freedom).  These 
effects  are  often  complex  in  nature  due  to  the  variable  and 
nonlinear  muscle  configurations  during  movement. 
Considering  only  the  effective  elastic  behavior  at  the  joint,  we 
can  regard  the  whole  muscle  elasticity  mapped  to  the  joint  space 
as  a  nonlinear  rotational  spring  attached  to,  and  “felt”  by  each 
joint.  Indeed,  there  exists  a  resultant  rotational  joint  stiffness  at 
each  joint,  which  has  the  same  effect  as  actual  muscle  elasticity. 
Therefore,  any  change  of  the  joint  angle  from  neutral  position 

will  result  in  a  restoring  torque  ^estonng  ?  which  can  be  linearly 
approximated  as: 

,  /'  =  (6) 

where  £7  is  appropriate  equivalent  rotational  joint  stiffness  for 
each  generalized  joint  spring,  and  qf  is  the  neutral  joint 

variable  corresponding  to  zth  joint  angle  qt .  The  neutral  joint 

variable  is  the  natural  unstretched  joint  configuration  where  the 
moment  about  a  joint  due  to  all  the  corresponding  muscle  forces 
are  balanced  without  any  external  load  or  joint  actuator  torque. 
The  coefficient  for  each  qt  is  obtained  from  the  simulation 

experiments  based  on  the  muscle  physiological  cross  sectional 
area  (PCSA).  In  vector-matrix  form,  this  equation  is  rewritten 
as  below. 


Restoring  =_K(q_q^  (7) 

where  is  the  vector  of  neutral  joint  angle,  and  K  is  a 
diagonal  matrix  of  joint  stiffness. 


LAGRANGE’S  EQUATIONS  OF  MOTION 

A  general  Lagrange’s  equation  in  vector-matrix  form  is 
d  d(W  +  T)  d(W  +  T)  _Q 
dt  dq  dq 

where  q  is  the  generalized  coordinate  vector,  T  is  the  total 
kinetic  energy  of  the  system,  W  is  the  virtual  work  done  by 
(non-inertial)  forces  and  torques.  The  virtual  work  is  the  sum  of 
the  virtual  works  done  by  conservative  forces  and  non¬ 
conservative  forces,  i.e.,  W  =  Wc  +  Wnc .  The  conservative  work 
can  be  expressed  as  the  negative  of  the  potential  energy  of  the 
force  system,  i.e.,  8Wc  =  -8V  .  Thus  we  can  write 

S(W  +  T)  =  ST  -  SV  +  SWnc  =  SL  +  SWnc  (9) 
where  V  is  the  total  potential  energy  of  the  system,  Wnc  is  the 
virtual  work  done  by  non-conservative  forces  and  torques,  and 
L  =  T -V  is  called  the  Lagrangian  function.  Then  equation  (8) 
can  be  written  as  follows. 

d_d_L  5L  d  3Wnc  8Wnc  Q 
dt  dq  dq  dt  dq  dq 


We  will  derive  more  explicit  form  of  the  equation  for  a 
general  serial  manipulator  and  force  system  shown  in  Figure  3. 

The  velocity  of  a  point  fixed  in  the  z'th  local  frame  with 
differential  mass  dm  in  link  z  can  be  expressed  in  terms  of 
global  coordinate  system  as  follows. 


dq. 


(11) 


J 


Then  the  total  kinetic  energy  K  of  the  human  link  system  is 


K  =  tJKi=\fj\Tr(yiy])dm  (12) 

**=1  Z  (=1  m 


Global  coordinate 


Figure  3.  Joint-link  system  with  external  loads 


3 


Copyright  ©  2005  by  ASME 


The  total  potential  energy  P  due  to  gravity  is  the  sum  of 
each  link’s  potential  energy. 

p  =  I^=I-^gr(Tn)  (13) 

7=1  7=1 

where  n  is  the  center  of  mass  vector  of  link  i  with  respect  to 
the  zth  local  coordinate  frame  and  g  is  the  augmented  (4x1) 
gravity  vector. 

Next,  we  assume  that  the  applied  non-conservative  forces 
are  actuating  torques  at  the  joints  and  external  loads  at  the  end- 
effector.  Then  the  variation  of  virtual  work  done  by  non¬ 
conservative  forces  is 


SWnc  =  Sq*T  +  Sx*F  +  £0. M  =  SqTr  - 


where  t  =  [r19r2 Tn]T  is  the  joint  torque  vector  (actuated 
by  muscles),  x  is  the  (3  x  1)  Cartesian  coordinate  of  the  point  of 

application  of  external  force  F  (a  3x1  vector),  and  0  is  the 
(3  x  1)  orientation  vector  of  the  link  with  respect  to  the  global 
frame,  where  the  external  moment  M  (a  3x1  vector)  is 
applied.  Assuming  that  the  velocity-dependent  forces  doesn’t 
exists  i.e.,  d(dWnc  / dq) / dt  =  0  ,  the  last  term  of  equation  (10) 
can  be  expanded  as 


where  J  is  the  Jacobian  matrix  (Sciavicco  and  Siciliano  1996). 
The  above  derivation  can  easily  extended  to  the  case  where 
multiple  external  loads  are  applied  to  any  location  of  any  link, 
not  necessarily  to  the  end-effector.  Likewise,  the  restoring 

torque  xRestormg  introduced  earlier  can  also  be  included  in  the 
above  equation  in  the  similar  manner. 

Based  on  the  procedures  given  in  Fu  et  al.  (1987)  with 
some  modifications,  equation  (10)  can  be  expanded  using  the 
kinetic  energy  (12),  the  potential  energy  (13),  and  the  extended 
form  of  non-conservative  work  (15).  Thus  the  final  vector- 
matrix  form  of  the  equations  of  motion  for  human  joint-link 
system  with  several  external  loads  is  shown  below  as  coupled, 
nonlinear,  second-order  ordinary  differential  equations. 


t  =  M(q)q  +  V(q,q)  +  £  Jf  m,g  +  £  J[ 


-K(q-qw)  (16) 


where  M(q)  is  the  mass-inertia  symmetric  matrix,  V(q,q)  is 
the  Coriolis  and  Centrifugal  vector,  IX  mi g  is  the  joint  torque 
vector  due  to  gravity  force,  J7  is  the  Jacobian  matrix  of  the 

position  vector  for  the  center  of  mass  of  the  /th  link,  and  Jk  is  the 
Jacobian  matrix  (for  any  point  in  any  local  coordinate)  of  the 


point  at 


(4x1)  location  vector  with  respect  to  kth  local 


coordinate  frame: 


~d°Tk(q)k  d°Tk(q)  k  d°Tk(q)  k  ] 

T  ,  \  n  L  •••  n  lk  •••  ~  lk  v1 ' ) 

Jfr(q)  “  d(h  d(lr  d(lk 

zo(q)  •••  Zz-i(q)  •••  z£-i(q)  _6 Xk 

where  Z._1?  i  =  is  the  local  z-axis  vector  of  joint  i 

expressed  in  terms  of  the  global  coordinate  system.  Here,  we 
only  take  the  first  three  elements  of  the  (4x1)  vectors 

(d°Tk(q)/dqj)krk . 

To  apply  the  above  equations  of  motion  to  the  human  joint- 
link  system,  some  physical  characteristics  of  anatomical 
segments  are  needed  to  determine  the  coefficients  in  the 
equations  of  motion.  These  physical  properties  are:  mass  of 
each  link,  center  of  mass  for  each  link,  moments  of  inertia  for 
each  link,  and  joint  stiffness. 

Our  digital  human,  Santos™,  is  modeled  based  on  a  typical 
male  soldier  anthropometry  data.  For  simplicity,  all  the  links  of 
Santos™  are  modeled  as  thin  rods,  so  that  the  moments  of 
inertia  are  calculated  easily  using  the  link  lengths  and  radii  of 
gyration.  The  mass  and  center  of  mass  data  are  obtained  from 
Chaffin  et  al  (1999)  and  Mi  (2004),  and  then  modified  to  fit 
into  our  Santos™  model.  The  joint  stiffness  is  estimated  from 
the  stiffness  of  all  the  muscles  involved  in  each  degree-of- 
freedom  joint  movement. 

HUMAN  ENERGY  EXPENDITURE  RATE 

Skeletal  muscles  may  be  thought  of  as  the  biochemical 
machines  with  chemical  energy  stored  in  ATP  going  into  the 
muscles  converted  to  mechanical  work  and  heat  energy 
(Sherwood  1997,  Umberger,  et  al.  2002).  In  other  words,  the 
total  metabolic  energy  expenditure  will  be  transformed  mainly 
into  the  sum  of  work  done  by  the  joint  torques,  heat  energy 
dissipation,  and  basal  metabolic  energy.  In  the  case  of  static 
loading,  where  the  mechanical  work  done  by  muscle  is  zero,  the 
muscle  energy  is  all  dissipated  as  heat. 

The  mechanical  power  is  defined  as  the  product  of  joint 
torque  and  joint  velocity.  The  total  mechanical  power  of  the 
system  w  is  the  sum  of  the  mechanical  power  for  all  the  joints. 

^=ZM=I M  (18) 

7=1  7=1 

Note  that,  since  the  torque  r.  and  the  joint  velocity  qi  are 

functions  of  time,  the  total  mechanical  power  W  is  also  a 
function  of  time.  In  other  words,  the  values  of  mechanical 
power  can  be  evaluated  at  each  time  instant.  The  absolute 
values  in  equation  (18)  mean  that  the  mechanical  energy 
consumed  by  the  muscle  cannot  be  reduced  by  simultaneous 
production  of  negative  power  by  another  muscle  at  the  same 
time,  or  by  the  same  muscle  during  different  periods  of  time. 

Based  on  several  heat  energy  formula  for  longitudinal 
muscle  model  (Hase  and  Yamazaki  1997,  Khang  and  Zajac 
1989),  the  muscle  heat  expenditure  is  divided  into  two 
components:  (i)  Muscle  maintenance  heat  and  (ii)  muscle 
shortening  heat.  The  muscle  maintenance  heat  is  the  energy 
consumed  in  proportion  to  the  muscular  tension.  The  muscle 
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shortening  heat  is  the  energy  consumed  in  proportion  to  muscle 
shortening  rate.  It  is  necessary  to  convert  these  relations  in 
muscle  space  into  the  generalized  coordinate  system,  i.e.,  the 
joint  space,  by  using  the  heat  energy  equation  given  in  Hase  and 
Yamazaki  (1997). 

For  a  given  muscle,  the  muscle  maintenance  heat  rate  is 
calculated  as  the  product  of  a  constant  and  muscle  activation. 
From  the  studies  on  muscle  activation  (Anderson  and  Pandy 
2001),  we  observe  that  the  generalized  maintenance  heat  rate  is 
approximately  proportional  to  the  joint  torque  and  inversely 
proportional  to  the  maximum  torque  limit.  So,  in  terms  of  joint 
space,  the  generalized  maintenance  heat  rate  Qm  is 


Qm  “IXK 

i= 1 


(19) 


where  hlm  is  the  generalized  coefficient  of  the 

maintenance  heat  at  joint  z.  These  coefficients  are  assigned  to 
each  joint,  and  are  obtained  from  the  maximum  voluntary 
torque  data. 

The  amount  of  muscle  shortening  heat  is  usually  small 
compared  to  that  of  the  maintenance  heat  (Bhargava,  et  al. 
2004,  Anderson  and  Pandy  2001).  For  this  reason  and  for 
further  simplicity,  we  neglect  the  contribution  of  muscle 
shortening  heat.  Furthermore,  this  term  resembles  the  form  of 
muscle  mechanical  power.  So,  by  adjusting  the  coefficients 
hlm  (z  =  1,  ...,z?)  of  the  generalized  maintenance  heat  Qm  ,  we 
can  obtain  good  approximation  of  total  muscle  energy  rate. 

The  basal  metabolic  rate,  BMR  is  defined  as  the  rate  at 
which  heat  is  produced  by  an  individual  in  a  resting  state.  BMR 
indicates  the  minimum  amount  of  energy  needed  to  carry  out 
metabolism  or  visceral  activity  at  rest.  We  use  the  formula  used 
by  Hase  and  Yamazaki  (1997)  for  estimation  of  BMR.  The 
BMR  of  the  whole  body  is  estimated  as  follows. 

B  =  0.685 W  +  29.8  (Joule/second  or  Watt)  (20) 

where  W  (Kg)  is  the  body  mass. 

Finally,  the  human  metabolic  energy  expenditure  rate  is  the 
sum  of  mechanical  power,  heat  rate,  and  BMR. 

^ Metabolic  *ZI  *)(0?/(0|  +  lX|r((0|  +  5  (21) 

7=1  7=1 


OPTIMIZATION  AND  MOTION/POSTURE  PREDICTION 

We  use  an  optimization-based  technique  for  motion 
prediction.  This  is  based  on  our  assumption  that  human  moves 
in  such  a  way  to  minimize  certain  human  performance 
measures,  subject  to  some  physical  and  physiological 
constraints.  Several  human  performance  measures  have  been 
investigated  and  shown  to  produce  various  natural  motions  and 
postures  (Yang  et  al.  2004,  Kim  et  al.  2004).  In  this  paper,  the 
metabolic  energy  expenditure,  as  described  in  the  previous 
section,  is  used  as  a  human  performance  measure  for  general 
motion  prediction. 

For  motion  prediction,  we  use  the  B- Spline  method,  which 
is  widely  used  for  manipulator  motion  studies.  Therefore,  the 


joint  angles,  joint  velocities,  and  joint  accelerations  are  all 
linear  functions  of  the  control  points  of  the  B-Spline  curves. 
The  optimization  procedure  calculates  the  joint  angle  profiles  of 
every  joint  in  the  form  of  B-Spline  curves  for  natural  upper- 
body  motion.  Also,  the  joint  velocity  and  acceleration  profiles 
are  directly  obtained  in  terms  of  control  points  from  the  direct 
calculation  of  the  B-Spline.  The  optimization  problem  for 
motion  prediction  can  be  stated  as  follows: 

Find:  Control  Points  ( P  z'  =  l,...,m;  z=l,...,z?) 

To  minimize:  Metabolic  Energy  ( EMetabolic ) 

Subject  to: 

Joint  limits  ( qf  <  q.  <  q1/  i  = 

Torque  limits  ( rf  <  ti  <  t]j  i  =  1, n ) 

Path  constraints  ( 1 1  x(q(zt))  -  p(/)  1 1  <  s ) 

Equations  of  motion 

T  =  M(q)q  +  V(q,  q)  +  £  J>,g  +  £  JtTF,  +  K  (q  -  qw  ) 

When  the  optimization  is  performed  at  each  time  interval, 
the  resulting  joint  kinematic  profiles  are  used  as  inputs  for  the 
equations  of  motion.  This  will  then  solve  the  inverse  dynamics 
problem  for  the  required  joint  torques.  The  energy  expenditure 
rate  at  each  time  instant  is  also  calculated  throughout  the 
motion. 

In  our  study,  we  use  SNOPT  software  package  developed 
by  Gill  et  al.  (2002).  The  SNOPT  package  is  a  general-purpose 
system  for  solving  optimization  problems  involving  many 
variables  and  constraints,  based  on  a  sequential  quadratic 
programming  (SQP)  method. 

EXAMPLE  RESULTS  AND  ANALYSIS 

In  this  section,  some  examples  of  lifting-motion  prediction 
and  joint  torque  evaluation  are  considered  for  the  21 -DOF 
Santos™  upper-body  model. 


Example  1 .  Prediction  and  analysis  of  dynamic  pulling  motions 
The  task  is  to  predict  the  motion  of  pulling  a  rope  with 
constant  tension  from  a  given  initial  position,  to  a  given  final 
position  in  2  seconds.  The  rope  is  stretched  horizontally  toward 
the  right  side  of  a  human  where  a  counter-weight  is  hanging  at 
the  other  end  of  the  rope  around  a  pulley.  Figures  4  and  8  show 
the  resulting  predicted  motions.  For  motion  1A  (Figure  4),  a 
0.001  Kg.  counter- weight  is  hanging  at  the  other  end  of  the 


rope. 

The  resulting  optimum  joint  kinematic  profiles  are  obtained 
at  each  0.04  seconds  and  are  shown  in  Figure  5.  At  each  time, 
the  joint  torque  values  and  the  metabolic  energy  rate  are 
calculated  (Figures  6  and  7).  The  metabolic  energy  expended 
for  this  motion  is  234.89  (Joules). 
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(a)  t  =  0  second  (initial)  (b)  t  =  2  seconds  (final) 
Figure  4.  Predicted  motion  of  Example  1A: 

Pulling  a  rope  with  0.001  Kg.  load 


Figure  5.  Predicted  joint  profiles  for  motion  Example  1A 


Figure  6.  Predicted  joint  torques  for  Example  1A 


Energy  Rate  Profiles  for  Example  1 A 


Figure  7.  Predicted  energy  rate  profile  for  motion  Example  1 A 


For  motion  IB  (Figure  8),  a  27  Kg  counter- weight  is 
hanging  at  the  other  end  of  the  rope.  This  task  is  same  as  1A, 
except  for  the  amount  of  counter- weight  (0.001  Kg  vs.  27  Kg). 
It  is  observed  that  Santos™  generated  different  motions  for 
different  amount  of  external  loads  required.  For  larger  tension 
(27  Kg.  weight,  in  this  case),  Santos™  tilted  his  trunk  to  use  his 
own  body  weight.  In  this  way  digital  human  could  also 
straighten  his  right  arm  to  minimize  the  joint  torques  at  the  wrist 
and  elbow. 


(a)  t  =  0  second  (initial)  (b)  t  =  2  seconds  (final) 

Figure  8.  Predicted  motion  of  Example  IB: 
Pulling  a  rope  with  27  Kg.  load 


Joint  Torque  Profiles  for  Example  IB 
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Figure  9.  Predicted  joint  torques  for  Example  IB 


In  Figure  9,  the  high  negative  torque  values  at  joints  1,  4,  7, 
and  10  for  motion  IB  indicate  the  major  contributions  of  these 
joints  to  the  left  lateral  bending  motion.  These  large  torques  are 
used  to  move  the  body  while  balancing  with  the  large  external 
pulling  force. 

The  metabolic  energy  consumed  during  this  motion  is 
614.24  (Joules).  As  expected,  this  value  is  much  higher  than 
the  energy  expenditure  of  pulling  with  0.01  N  (234.89  (Joules)). 


Example  2.  Lifting  same  weight  at  different  locations 

In  this  and  the  following  examples,  we  consider  the  initial 
posture  of  each  motion  to  predict  the  risk  level  of  injury  at  the 
onset  of  the  task. 

Figures  10  and  11  show  two  different  onset  postures  of 
manual  lifting  tasks.  In  posture  2A,  a  9  kg  object  on  hand  is 
located  far  right  from  the  body  and  in  posture  2B,  the  same 
weight  is  located  far  front  from  the  body. 
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With  the  velocities  and  accelerations  set  to  zero  (static 
onset  posture),  the  equations  of  motion  (16)  are  used  to  obtain 
the  torque  at  each  joint.  The  joint  torques  of  both  postures  are 
plotted  in  Figure  12. 


Figure  10.  Posture  2 A  Figure  11.  Posture  2B 


Joint  Torques  of  Lifting  Postures 


100 

50 
0 

-50 
-100 

-150 

Joint  Number 


Figure  12.  Joint  torques  of  lifting  postures  2 A  and  2B 

In  posture  2 A,  the  spine  joints  for  lateral  bending  exert 
large  torques  for  right  lateral  bending,  while  in  posture  2B,  the 
spine  joints  for  flexion/extension  exert  large  torques  for  flexion. 
The  large  torque  values  for  different  motion  segments  will  have 
different  effects  on  the  body  and  may  cause  different  kind  and 
risk  level  of  injury  (especially  on  the  low-back).  It  can  be 
observed  that  the  two  postures  have  different  torque 
distributions  at  the  shoulder  and  arm  joints. 

Using  equation  (21),  the  energy  expenditure  rate  is 
calculated  as  275.5  (Joules/sec)  for  posture  2A  and  308.3 
(Joules/sec)  for  posture  2B.  These  values  are  within  reasonable 
range  of  energy  expenditure  rate  during  weight  lifting  exercise 
(Sherwood  1997). 

Example  3 .  Lifting  same  weight  with  different  postures 

A  manual  task  of  lifting  a  10-Kg  object  with  two  different 
postures  is  studied  in  this  example.  The  two  predicted  postures 
are  shown  in  Figures  13  and  14.  In  each  case,  Santos™  is 
holding  the  object  so  that  the  palm  faces  upward.  Posture  3 A  is 
the  lifting  posture  where  the  right  arm  is  stretched  forward  and 
the  object  is  far  from  the  body,  located  at  (-30,  20,  65)  (cm) 
with  respect  to  the  global  frame  at  the  lumbar  region  (Figure  1). 
In  posture  3B,  the  right  elbow  is  bent  and  the  object  is  brought 
closer  to  the  body,  i.e.,  (-30,  15,  40)  (cm)  from  the  global  frame. 


□  Posture  2A 
■  Posture2B 


The  resulting  joint  torques  for  both  postures  are  plotted  in 
Figure  15. 


(a)  Oblique  view  (b)  Side  view 
Figure  13.  Posture  3  A 


(a)  Oblique  view  (b)  Side  view 
Figure  14.  Posture  3B 


Joint  Torques  of  Lifting  Postures 


Joint  Number 


□  Posture3A 
■  Posture3B 


Figure  15.  Joint  torques  of  lifting  postures  3  A  and  3B 


As  indicated  in  Figure  15,  holding  the  object  far  from  the 
body  requires  the  muscles  of  the  low  back  to  exert  an  extension 
torque  of  71.7  (N.m)  at  joint  2  (flexion/ extension  joint  of  low- 
back).  However,  if  the  weight  is  located  closer  to  the  body, 
only  37.7  (N.m)  torque  is  required  at  joint  2.  The  similar 
change  can  be  observed  at  other  joints  that  exert  extension 
torques,  i.e.,  joints  5,  8,  and  11.  The  reduced  joint  torque 
exerted  by  the  back  muscles  will  diminish  the  load  experienced 
by  the  spine  and  therefore  will  reduce  the  risk  of  back  injury. 

The  big  difference  in  torque  values  can  be  found  also  in 
joints  14  and  16.  These  are  abduction/adduction  of  the  clavicle 
and  shoulder  joint,  respectively.  Just  as  the  low  back  joint  case, 
this  is  due  to  the  difference  in  the  lever-arm  length.  Thus, 
posture  3A  is  more  likely  to  generate  muscle  fatigue  than 
posture  3B.  This  gives  another  reason  that  the  load  should  be 
kept  close  to  the  body  during  manual  lifting. 


7 


Copyright  ©  2005  by  ASME 


Holding  an  object  at  a  far  distance  from  the  body  requires  a 
large  amount  of  torques  to  be  exerted  by  the  back  and  shoulder 
muscles  while  the  amount  of  torques  can  be  reduced  by  holding 
the  object  closer  to  the  body.  Thus,  our  result  supports  the 
recommendation  that  workers  should  keep  the  load  close  during 
manual  lifting  tasks  to  prevent  injury. 

The  energy  expenditure  rate  for  posture  3A  is  calculated  as 
288.3  (Joules/sec),  while  for  posture  3B,  it  is  230.4 
(Joules/sec).  It  can  be  seen  that  the  value  of  energy  expenditure 
rate  for  any  time  instant  gives  us  a  single  criteria  that  evaluates 
the  level  of  “effort”  of  a  given  posture. 

CONCLUSIONS  AND  DISCUSSION 

An  optimization-based  method  and  examples  of 
motion/posture  prediction  and  analysis  for  manual  lifting  were 
introduced.  The  results  have  shown  that  our  model, 
formulation,  and  method  can  be  extended  to  analyze  any  kind  of 
general  material  handling.  Our  method  for  motion  prediction 
and  analysis  can  be  very  helpful  in  determining  the  risk  of 
injury  during  manual  lifting  tasks.  The  evaluation  of  joint 
torques  during  lifting  can  help  the  user  to  redesign  tasks  so  that 
the  stresses  on  the  low-back  and  shoulder  are  decreased.  The 
simulation  also  showed  how  Santos™  generates  different 
motions  and  postures  in  response  to  different  amount  of 
external  loads. 

A  general  form  of  dynamic  equations  of  motion  were 
derived  and  used  to  obtain  the  metabolic  energy  expenditure 
formula.  It  is  also  used  to  solve  the  inverse  dynamics  problem 
of  the  motion  prediction.  This  equations  of  motion  include  the 
general  external  loads  (forces  and  moments)  applied  not  only  to 
the  end-effector  (hand),  but  also  to  any  other  parts  of  the  body. 
This  generalization  of  point  of  force-application  enables  us  to 
simulate  the  cases  when  a  human  is  hanging  a  bag  on  the 
forearm,  wearing  a  backpack  on  the  shoulder,  and  so  forth.  The 
metabolic  energy  expenditure  was  used  as  a  cost  function  to  be 
minimized.  By  rearranging  the  metabolic  energy  rate  formula, 
it  is  observed  that  the  energy  rate  is  a  summation  of  weighted 
joint  torques,  where  the  weight  function  is  the  sum  of  the 
corresponding  joint  velocity  and  the  heat  coefficient.  This 
minimum  energy  cost  criteria  implicitly  indicates  minimum 
joint  torques  and  thus,  the  motion  based  on  minimum  energy  is 
equivalent  to  the  motion  of  minimum  joint  torques. 

The  joint  torque  solutions  of  the  inverse  dynamics  can  be 
used  as  an  input  for  the  prediction  of  muscle  force  distribution. 
For  the  purpose  of  motion/posture  prediction,  the  principle  of 
generalized  coordinates  and  generalized  torques  provided  us 
with  a  simplified  and  accurate  kinematical  model.  However, 
multiple  muscles  contribute  to  a  single  degree-of-freedom  joint 
motion  and  the  actual  configuration  of  muscles  varies  as  the 
corresponding  joint  moves.  This  leads  us  to  the  future  studies 
of  detailed  biomechanical  analysis  on  muscle  forces  and 
deformations,  and  then  to  muscle  fatigue  prediction.  Finally,  to 
obtain  more  accurate  value  for  the  energy  expenditure  and  joint 
torques,  we  will  acquire  more  exact  data  of  physical  properties, 


such  as  centers  of  masses,  moments  of  inertia,  and  products  of 
inertia. 

Although  our  simulation  results  are  kinematically  realistic, 
the  actual  human  body  is  more  complex  and  flexible  than  the 
Santos™  model.  Especially  for  the  detailed  study  of  the  spine, 
a  higher  degrees-of-freedom  model  should  be  used  to  improve 
accuracy.  Our  future  work  is  to  incorporate  muscles  to  build 
higher  degrees-of-freedom  model  that  is  anatomically  realistic. 
This  will  also  enable  us  to  investigate  the  muscle  forces  and 
fatigue  properties  during  material  handling  task. 
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